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Abstract 

We derive semiclassical expressions for the Kubo conductivity tensor. Within 
our approach the oscillatory parts of the diagonal and Hall conductivity are 
given as sums over contributions from classical periodic orbits in close relation 
to Gutzwiller's trace formula for the density of states. Taking into account 
the effects of weak disorder and temperature we reproduce recently observed 
anomalous phase coherence oscillations in the conductivity of large antidot 
arrays. 

03.65Sq, 73.50.Jt 



Typeset using REVTeX 



* Unite de Recherche des Universites Paris XI et Paris VI associee au CNRS 



1 



Interference phenomena in quantum electron transport through small microstructures 
are usually interpreted within two complementary frameworks: The Landauer-Biittiker for- 
malism commonly describes transmission through single phase coherent devices 0. The 
current is expressed in terms of (a sum over) conductance coefficients between channels in 
different leads attached. On the other hand linear response theory (Kubo formalism) has 
proved useful to treat bulk transport properties of samples with a size exceeding the phase 
breaking length. Besides coherence effects related to the presence of disorder (e.g., universal 
conductance fluctuations, weak localization) the development of high mobility devices has 
opened experimental access to the ballistic limit where the elastic mean free path is consid- 
erably larger than the system size and the conductance reflects the geometry or potential 
landscape of the microstructures. This has especially oriented interest to questions how clas- 
sically regular or chaotic electron dynamics manifests itself on the level of quantum transport 
0-^]. In this spirit a semiclassical approach to conductance within the Landauer-Biittiker 
framework has already been successfully performed expressing the conductance coefficients 
in terms of interfering electron paths [[|. 

Recent experiments on magnetotransport in antidot structures unveiled the lack of a cor- 
responding semiclassical approach to Kubo bulk conductivity Antidot superlattices 
consist of arrays of periodically arranged holes " drilled" through a two dimensional electron 
gas (2DEG). Since the lattice constants a (of typically 200-300 nm) are significantly larger 
than the Fermi wavelength ~ 50 nm the dynamics of the electrons moving in between 
the repulsive antidots can be considered to be of semiclassical nature. The combined poten- 
tial of the superlattice and a perpendicular magnetic field gives rise to a variety of peculiar 
effects: The diagonal magnetoresistivity p xx exhibits pronounced peaks due to trapping of 
electrons encompassing a particular number of antidots at magnetic field strengths related 
to specifc values of the normalized cyclotron radius R c /a J7|. Superimposed upon these 
resistivity peaks (reflecting the classical chaotic electron dynamics of the antidot geometry 
||) additional quantum resistivity oscillations of anomalous periodicity had been observed 
at sufficiently low temperature indicating phase coherence phenomena PH. They cannot 
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be attributed to interfering electron waves traversing the whole device (as in the case of 
single junctions) since the antidot arrays are too large to maintain phase coherence. How- 
ever, assuming that p xx reflects density of state oscillations we were able to describe the 
periodicity of the modulations observed in terms of quantized periodic orbits in the antidot 
array Nevertheless, a complete direct semiclassical approach to the conductivity tensor 
was still missing || in spite of important work by Wilkinson providing a semiclassical eval- 
uation of the diagonal conductivity in a somewhat different context ||10|| . Such an approach 
is also desirable since the antidot measurements have not yet been completely reproduced 
by quantum mechanical calculations which turn out to be rather involved fTT|j . 

In this letter we derive semiclassical Kubo-type transport formulas by stationary phase 
evaluation of the disorder averaged two-particle Green functions. We obtain diagonal and 
Hall conductivity oscillations in terms of sums over periodic orbit contributions with the 
classical actions of the orbits determining the periodicities. Using a model antidot potential 
and working at finite temperature we are able to reproduce the amplitudes and periodicities 
of measured quantum oscillations in antidot superlattices. 

Within the Kubo formalism the static conductivity tensor at temperature T is given by 



(1) 



where V is the total area, g s accounts for the spin degeneracy, f(E) = 1/[1+ exp((E— E F )f3)] 
is the Fermi function (/3 = 1/kT) and <> denotes an average over weak disorder [jnj. The 
diagonal and Hall conductivities = <Jij(E) can be written as [T3| 



a xx = 7re 2 ftTr {vJ(E - H)vJ(E - H)} (2a) 
<T*v = e 9N( f B B) + ^ Tr {v*G + (E)v y 5(E -H)- vJ(E - H)v y G~(E)} . (2b) 

Here, the Vi are the operators of the velocity components, G + (E) (G~(E)) is the advanced 
(retarded) Green function, N(E) = J e Tt5(E' — H)dE' the number of states below E and 
B is the magnetic field. 

For the semiclassical evaluation of Eq. (0) it is convenient to use 5{E — H) = — [G + (E) — 



G (E)]/2m and to work within position representation. Semiclassically, the operators V{ 
act as classical quantities and with Gg(r, r') = G^*(r', r) we can reduce Eq. (Q) to 

0xx = ~2ir Re {J dr J dr'v x G%(r,r')v' x G%*(r, r') - G%(r', r) | , (3a) 
^ = e^P + ^Re{| dr | dr'u^r, r')< [G£*(r, r') - G+(r', r)]} . (3b) 

We write the Green function as the Laplace transform of the Feynman propagator which is 
semiclassically given in terms of a sum over all classical paths a from r' at time t' — to r 
at time t @: 

% 7T 

-W a (r',0;r,t) + m a i- 
n 2 



G+(r,r') = -—j / die***/* ]T A,(r', r) exp 



Tin Jo a Lft 2. 

Wo, = JgCdt is the Hamilton principal function D a (r', r) = |detjj(— dp^^v', 0; r, t)/^)) 1 / 2 
contains information of the phase space in the vicinity of the trajectory a and m a is a 
topological phase. The factor exp(— t/2r s ) (r s being the elastic scattering time) results 
from calculating the disorder averaged Green functions within Born approximation ||15|| . 
Inserting the representation Eq. ([|) of G% and G^* into Eq. (^) leads to double sums J2 a & 
of products of path contributions. We begin with the semiclassical evaluation of the traces 
over the terms ViGg(r, r')v'jG% (r, r') in Eq. (^). The phases of the diagonal parts (a = a) 
cancel out and we obtain the classical contribution to the conductivity in two dimensions 
after one time integration and a transformation of the integral over the final coordinate 
into an integral over initial momenta (by means of D a which acts as a Jacobian), i.e., 
a ij = e 2 /h 2 J Q °° dt < Vi(t)vj(Q) > pq exp(— t/r 3 ) . Here, <> pq denotes a classical phase space 
average over the energy shell. Such a Kubo formula had been used in Ref. || to calculate 
numerically the classical part of the conductivity for antidot arrays. 

In order to take care of quantum phase coherence phenomena we have to take into account 
non-diagonal (a ^ a) parts in the sum J2 a a- Consider, e.g., the contribution of Green 
function products of direct paths a (from r' to r) and paths a with at least one selfcrossing 
as depicted in Fig. 0(a). We will compute the trace integrals — as usual for a semiclassical 
approach — by stationary phase approximation. The stationary phase condition for the r 



integral requires for the final momenta p£ = from which follows that the paths a and a 
must coincide between r' and r (Fig. 0(b)). Stationary phase approximation for the trace 
integral over r' selects all pairs of paths with the same initial momenta, i.e. = p^. Both 
conditions are fulfilled for all pairs of orbits with initial and final phase space points lying on 
a periodic orbit (Fig. |l](c)) . Besides the example shown in Fig. |l](c) we must include all pairs 
of orbits with the shorter one showing an arbitrary number of full traversals before reaching 
the final point r and the longer one having r additional repetitions along the periodic orbit. 
The semiclassical evaluation of the trace integrals and one time integral follows similar 
lines as Gutzwiller's derivation of the periodic orbit trace formula for the density of states 
I^ , P^6|JT7| . After performing finally the convolution (Eq. ([[])) with the Fermi function we 



obtain from the ViGg v'a G%* terms oscillating parts of the diagonal and Hall conductivity as 

sums over contributions from all primitive periodic orbits (po) of the classical system and 

their repetitions r: 

2a e 1 00 R (3) P ^T^/(2r a ) 

a° s x c (E F , B- T) = V CZ V m cos [r(S po /h - //V/2)] , 

xxK ' V h % xx f^ x |det(Mr-l)| 1/2 

In Eq. (|) S P °(E,B) = § po Y>d<\ is the classical action, T po the period and /i po the Morse 
index of a periodic orbit. The monodromy matrix Nl po measures its classical stability (16 



/■oo , r TP° 

Cf° = J dte~ t/Ts J o dTVi(r)vj(t + T) (6) 

are velocity-correlation functions for motion along the periodic orbits. The function 

; t c = — (7) 

smh(ri t °/Tp) 71 

gives rise to temperature damping which is exponential for orbits with T po ^> Tp. 

Stationary phase evaluation of the trace integrals in Eq. (fj) for the terms 
ViGg(r, r')u'-G^(r', r) representing recurring paths from r via r' leads also to periodic or- 
bits (see Fig. |l](d-f)). However, the corresponding velocity correlation functions vanish and 
therefore ViG + v'jG + terms do not contribute to semiclassically. 



Up to the correlation functions C v x ° x the trace formula Eq. (|5a| ) for a xx is exactly the 
same as Gutzwiller's trace formula for the density of states, quantum oscillations are 
semiclassically related to interference (due to the phase differences ~ rS po /h) between pairs 
of paths along a periodic orbit differing by r in their number of full traversals. This approach 
therefore requires only phase coherence on the length scale of a periodic orbit and not 
throughout the entire system. 

In the following we will apply our results to magnetotransport in 2DEG's. In the unmod- 
ulated case the electrons follow cyclotron orbits with action S cyc = 2ttEf/uj and frequency 
to = eB/m*. Then the classical amplitudes 0] and velocity correlation functions in the 
trace formulas ([5]) can be calculated analytically. For the diagonal conductivity we obtain, 
e.g., 

(n s being the carrier density) which represents a semiclassical approximation of the Shub- 
nikov de Haas oscillations and coincides with the quantum mechanical result for a constant 
scattering time t s [I5|. In the case of a high mobility 2DEG (ur s ^> 1) with a distinct 



(8) 



Landau level structure the use of a density of state dependent scattering time is usually 
required in order to obtain the correct conductivity amplitudes. However, in the antidot 
arrays which we consider in the following the periodic superlattice potential strongly mixes 
the Landau levels [llj and therefore justifies the use of an energy independent r s . The inset 
of Fig. H|(a) depicts the measured magnetoresistivity p xx (B) of an antidot array for two tem- 
peratures T=4.7 K and 0.4 K. While the gross structure of the broad (classical) maximum 
at 2R C (B) = a persists up to temperatures of 60 K the (quantum) oscillations at T=0.4 K 
on top of the maximum disappear with increasing temperature (at T=4.7 K). 

In order to show that the quantum oscillations at T=0.4 K are related to a x s x (B) we 
compare the experimental difference signal Aa xx = a xx (T=0AK)—a xx (T=4.7K) shown in 
Fig. g(a) with our semiclassical predictions for a° s x (B). As a model of the antidot lattice 
we use the potential V(x, y) = V \ sin(7ra;/a) sin(7n//a)| /3 which has already proved useful to 



describe magnetotransport in antidot arrays |3]||. Vq is determined through the normalized 
width d/a of the antidots at the Fermi energy and f3 (governing the steepness of the antidots) 
remains as the only free parameter. Here we use d/a = 0.5 and (3 = 2.3. 

In order to calculate numerically a° s x (B) it is convenient |20] to expand the prefactors 



in Eq. Q5*a| ) into geometrical series and to perform the sum over r to get the following form 
a~(B r , B; T) = f S - Re £ Cgi- £ — «W (9) 



with 



t M = (±i)* exp 



i ( x -^)"(* + 5 A- - (2/ + i; 



(10) 



T/3 2r s 

Here, A po = i^ po ) 7 po (real) being the winding number, k a semiclassical quantum number 
for a stable periodic orbit and X po > (real) being the Liapunov exponent in case of an 
unstable periodic orbit |2(J. The "— " sign in Eq. ( |10"D applies to unstable inverse hyperbolic 



orbits; otherwise the "+" sign holds. A correct application of Eq. @ requires a detailed 
knowledge of the phase space structure of the underlying classical motion in the potential of 
the antidots and the magnetic field. It turns out that the major part of the classical phase 
space is chaotic with a few islands of stable motion imbedded which vanish with decreasing 



magnetic field [^l]]. Under the experimental conditions (T = 0.4 K and ujt s pa 2 at 2R C = a 
f2~2fl ) only the shortest periodic orbits contribute significantly to a° s x c since the terms from 
longer orbits (T po > r s or T po > Tp) are exponentially small. For our calculations we take 
into account the seven fundamental periodic orbits shown in Fig. |2|(d). 

Our result for a° s x c (at T=0.4 K) is shown as the full line in Fig. @(b): a° s x c (B) oscillates 
with the same frequency as the measured conductivity. The period of the oscillations is 
nearly constant with respect to B in contrast to the 1/B periodic behaviour of ordinary 
Shubnikov de Haas oscillations (Eq. (§)) which are shown in Fig. ^|(c) for an unmodulated 
2DEG under the same conditions. The B periodic behaviour of quantum oscillations in 
certain antidot arrays can be related to the mechanism that flux enclosing periodic orbits 
between antidots cannot expand with decreasing B as in the case of free cyclotron motion 



0. Having in mind that the antidot potentials are not precisely known also the semiclas- 
sical amplitudes are on the whole in agreement with the measured curve of Aa xx . They 
show — as in the experiment — an irregular behaviour resulting from interference effects 
between different periodic orbits and from the magnetic field dependence of the classical or- 
bit parameter entering Eq. The semiclassical curves show also the correct temperature 
dependence: Due to the temperature damping R T (f3) (Eq. (|7|)) the periodic orbit oscillations 
decrease for T = 2.5 K (dotted line in Fig. §(b)) and nearly disappear at T = 4.7 K (dashed 
curve) as in the experiment (dashed curve in the inset of Fig. 0(a)). Fig. 0(e) illustrates the 
ojt s dependence of o° xx (B) for three values ujt s — 1, 2 and 5 showing that the result does 
not strongly depend on the r s chosen. 

In summary we derived semiclassical formulas for the magnetoconductivity tensor within 
the Kubo approach and applied them to quantum transport in antidot lattices. However, 
it remains to understand the semiclassical relation between the Kubo formalism leading 
to periodic orbit contributions and the semiclassical Landauer Biittiker approach based on 
interference effects due to non-closed classical trajectories since the two corresponding quan- 
tum mechanical approaches are known to be equivalent under certain conditions 
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FIGURES 

FIG. 1. Typical semiclassical paths a, a contributing to (a) ViGg(r,r')vjGg*(r,r') and (d) 
ViGfl(r, r')v'jGfl(r' , r). The stationary phase condition at r selects paths with p Q (r) = p«(r) (b,e). 
Paths along periodic orbits result from the stationary phase condition p a (r') = Po(r') at r' (c,f). 

FIG. 2. (a) Oscillatory part Aa xx = a xx (0AK) — a xx {^.1K) of the experimental M diagonal 
conductivity (a xx = p xx /(p xx + P xy )) of an antidot array as a function of an applied magnetic 
field. Inset: Measured resistivity p xx at T = 0.4 K (full line) and T = 4.7 K (dashed line), (b) 
Semiclassically calculated oscillatory part a° s x (from Eq. |||) for three different temperatures T = 
0.4 K (full line), T = 2.5 K (dotted line) and T = 4.7 K (dashed line), (c) Semiclassical Shubnikov 
de Haas oscillations (Eq. |8|) for an unmodulated 2DEG under the same conditions as in (b). (d) 
fundamental periodic orbits in a model antidot potential which enter the semiclassical calculation, 
(e) a° x s x c at T = 0.4 K for ojt s = 5 (dotted line), 2 (full line) and 1 (dashed line). 
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